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We investigate the possibility that dark matter is a mixture of two non-interacting perfect fluids, 
with different four-velocities and thermodynamic parameters. The two-fluid model can be described 
as an effective single anisotropic fluid, with distinct radial and tangential pressures. The basic equa- 
tions describing the equilibrium structure of the two-fluid dark matter model, and of the tangential 
velocity of test particles in stable circular orbits, are obtained for the case of a spherically symmet- 
ric static geometry. By assuming a non-relativistic kinetic model for the dark matter particles, the 
density profile and the tangential velocity of the dark matter mixture are obtained by numerically 
integrating the gravitational field equations. The cosmological implications of the model are also 
briefly considered, and it is shown that the anisotropic two-fluid model isotropizes in the large time 
limit. 



PACS numbers: 04.50. Kd, 04.20. Cv, 04.20. Fy 



> 
(N 

(N 

o 



X 



INTRODUCTION 



Cosmological observations provide compelling evi- 
dence that about 95% of the content of the Universe re- 
sides in two unknown forms, denoted dark matter and 
dark energy, respectively. The former resides in bound 
objects as non-luminous matter, and the latter in the 
form of a zero-point energy that pervades the whole Uni- 
verse [l|. Dark matter is thought to be composed of 
cold neutral weakly interacting massive particles, beyond 
those existing in the Standard Model of Particle Physics, 
and not yet detected in accelerators or in direct and in- 
direct searches. There are many possible candidates for 
dark matter, the most popular ones being the axions and 
the weakly interacting massive particles (WIMP) (for a 
review of the particle physics aspects of dark matter see 
) . Their interaction cross section with normal baryonic 
matter, while extremely small, are expected to be non- 
zero and we may expect to detect them directly. Scalar 
fields or other long range coherent fields coupled to grav- 
ity have also been used intensively to model galactic dark 
matter Q ■ Recently, the possibility that the galactic dy- 
namics of massive test particles may be understood in 
the context of modified theories of gravity, without the 
need for dark matter, was also explored [4]. 

Despite these important achievements, at galactic 
scales ~ 10 kpc, the ACDM model meets with severe 
difficulties in explaining the observed distribution of the 
invisible matter around the luminous one. In order for 
galactic gas and stars to attain a constant velocity in 
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the outer regions, the density profile of the dark mat- 
ter must fall approximately with r~ 2 in the respective 
regions. These profiles resemble the density structure 
of an isothermal and self-gravitating system of particles, 
characterized by a constant velocity dispersion. In fact, 
A^-body simulations, performed in this scenario, predict 
that bound halos surrounding galaxies must have very 
characteristic densityprofiles that feature a well pro- 
nounced central cusp [5J . On the observational side, high- 
resolution rotation curves show, instead, that the actual 
distribution of dark matter is much shallower than the 
above, and it presents a nearly constant density core @. 

The possibility that dark matter could be described by 
a fluid with non-zero effective pressure was also consid- 
ered in the physical literature Q. In particular, in Q it 
was assumed that the equation of state of the dark mat- 
ter halos is polytropic. The fit with a polytropic dark 
halo improves the velocity dispersion profiles. 

Due to their importance in astrophysics and cosmol- 
ogy, two-fluid models have been extensively investigated 
in the physical literature. In particular, dust shells of fi- 
nite size, moving with different four-velocities and form- 
ing a common region, where the two dust components ex- 
ist simultaneously, were considered in ll| for the case of a 
Lemaitre-Tolman geometry. To model voids a spherically 
symmetric dust with different energy density values has 
been matched by an intermediate region. By introducing 
two dust components the boundaries no longer need to 
be comoving. The two-component dust model was gen- 
eralized in [121 ] by introducing an additional homothctic 
vector. The obtained metrics are regular, except for a 
big bang singularity and in the limit of large times they 
can be given by an asymptotic expansion of the metric 
coefficients. Using this family of solutions, voids whose 
edge is not comoving in comparison with its surroundings 
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can be modeled in geometric terms. 

In this work, we consider the possibility that dark mat- 
ter may be modeled as a mixture of two non-interacting 
perfect fluids, with different four-velocities. This con- 
figuration is formally equivalent to a single anisotropic 
fluid [H, We obtain the basic equations describing 
the structure of this system, namely, the mass continuity 
and the hydrostatic equilibrium equations, as well as the 
equation giving the tangential velocity of test particles 
in stable circular orbits. For the dark matter we adopt 
a non-relativistic kinetic description, with the fluid pres- 
sure proportional to the density of the fluid and to the ve- 
locity dispersion of the constituent particles. By assum- 
ing that the two fluids have the same density, the dark 
matter profile and the tangential velocity are obtained by 
numerically integrating the field equations. Some cos- 
mological implications of the model are also discussed, 
and it is shown that in the large time limit the initially 
anisotropic homogeneous two-fluid dark matter distribu- 
tion ends in an isotropic phase. 



by an energy-momentum tensor having the form given by 
Eq. ([1} can be significantly simplified if we cast it into 
the standard form of perfect anisotropic fluids. This can 
be done by means of the transformations 0, [H| 




P2+P2 
P1+P1 



Hasina, (2) 



-> W*» = Mucosa- J Pl +Pl L^sina, (3) 
V 92 + P2 



representing a "rotation" of the velocity four-vectors in 
the (U^jW^) velocity space, which leave the quadratic 
form (pi + Pi) U^U v + {p2 + P 2 ) WW invariant. Thus, 
T» v (U,W) = T^ v (U*,W*). Next, we choose U*^ and 
W*^ such that one becomes timelike, while the other is 
spacelike. Therefore the two vectors satisfy the condition 
U*^W* = 0. With the use of Eqs. ©-© and U*»W* = 
we obtain the rotation angle as 



II. DARK MATTER AS A MIXTURE OF TWO 
PERFECT FLUIDS 



tan*. = 2^ {P1+Pl \ {p2+P jww,. (4) 

Pl+Pl — \P2 +P2) 



We start our study of dark matter by assuming that 
it consists of a mixture of two perfect fluids, with energy 
densities and pressures pi, Pi and P2, P2, respectively, 
and with four velocities C/ M and W^, respectively. The 
fluid is described by the total energy-momentum tensor 
T* 1 ", given by 



By defining the quantities V = U*^ / y/U* a U*, 

e = T^V^ = 

* = T^XvX* = 
and II = Pi + P 2 , re- 



X 



= W^/yZ-W^W*, 

( P i+Pi)U* a U* - (P1+P2 
(Pi + P2) - (P2 + P2) W* a W c 
spectively, the energy-momentum tensor of the two non 
interacting perfect fluids can be written as 



= (pi + Pi)EW -Pig' 1 " + 
( P2 + P 2 )W^W V -P 2 g^. 



(1) 



The four-velocities are normalized according to 
U^Uf, = I and W^W,! = 1, respectively. In the 
present paper we use the natural system of units with 
c = G = 1. The study of the physical systems described 



T» v = (e + H) VV - Tig* 1 " + - II) x"x" , (5) 

where V% = I = -x^Xp, and = HQ]}. The 

energy-momentum tensor given by Eq. §5§ is the standard 
form for anisotropic fluids fioj . 

The energy density e and the radial pressure \& are 
given by 



J 



"(Pi 



P2-P1-P2 



^\/(pi+Pl+ P2+P2) 2 + 4!/, M -Pj)I,m - P.) 



(6) 



{pi + P2-P1-P2 



Pi 



P2 



P2) 2 + 4 (pi + Pi) (p 2 + P 2 ) (UnWrf, 



(7) 



respectively (9I. Hfjj. 

In comoving spherical coordinates x° 



t, x 1 



— i9, and x 
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xx f may choose V 1 = V 2 = V 3 



r, 
0, 



V°V = I, and x° = X 2 = X 3 = 0, X X Xi = -iSES 



X = X 

Therefore the components of the energy-momentum of 
two non-interacting perfect fluids take the form: Tg = e, 
and Tf = T| = -II, where e is the total 



n 



energy- density of the mixture of fluids, ^ = P r is the 
pressure along the radial direction, while II = P± is the 
tangential pressure on the r = constant surface. 
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III. ANISOTROPIC FLUIDS IN SPHERICALLY 
SYMMETRIC STATIC SPACETIMES 



IV. DARK MATTER AS A MIXTURE OF TWO 
FLUIDS 



In the following we restrict our study of the two- 
component dark matter to the static and spherically sym- 
metric case, with the metric represented as 

ds 2 = e^dt 2 - e x ^dr 2 - r 2 (dd 2 + sin 2 M<p 2 ) . (8) 

For the metric given by Eq. ([8]). the Einstein gravi- 
tational field equations, describing the mixture of two 
fluids, can be easily integrated to give 

* _ . 2M(r) 



dv 
dr 



with M{r) = 4?r / er z dr, (9) 



AIL- * 



r e 



9 ' 



dM 
dr 



Aner 



(10) 



(11) 



d9 _ (e + (47r*r 3 + M) 2 
~dr~ ~ r 2 (1 - 2M/r) + r 



(n-*). (12) 



Equation (TTO]) is the consequence of the conservation of 
the energy-momentum tensor, T£ M = 0, while Eqs. (fTTj) 
and (jT2J) are the mass continuity, and the hydrostatic 
equilibrium (TOV) equation. 

The galactic rotation curves provide the most direct 
method of analyzing the gravitational field inside a spiral 
galaxy. The velocity of the cloud, as measured by an 
inertial observer far from the source, is given by 



dCl 
~dJ 



(13) 



For a stable circular orbit dr/dt = 0, and the tan- 
gential velocity of the test particles can be expressed as 
2 where dfl 2 = dd 2 + sin 2 iW0 2 . In 



J tg 



V (dfl/dt) 



terms of the conserved angular momentum I = r </) and 



energy E 



t/2, by v% 



e v ( r H, the tangential velocity is given, for 
e v l 2 jr 2 E 2 . By taking into account the 
explicit expressions for I and E we obtain for the tangen- 
tial velocity of a test particle in a stable circular orbit 
the expression v 2 g = ri/' /2. 

The tangential velocity profile can be obtained as a 
function of the total density and of the mass of the dark 
matter from Eq. (|TU| as 



vl B {r) 



4tdIt 3 + M 
r (1 - 2M/r) ' 



(14) 



Eq. (|12|) . can be written in the following equivalent 
form 



d® 
dr 



9) 



v 2 g (r) + -(LI-9) 



(15) 



The kinetic energy-momentum tensor T{? associated to 
the frozen distribution of dark matter is given by T{? — 
(d/h 3 ) J d 3 pf(p)p^p u /p° , where f(p) is the dark matter 
particle distribution function, p^ is the four-momentum, 
p is the three-momentum, with absolute value p, and g 
is the number of helicity states, respectively. The energy 
density e of the system is defined as 



(16) 



while the pressure of a system with an isotropic distribu- 
tion of momenta is given by 



P=±f I n,M* P = JLf£f< p )* Pi 



(17) 



where the velocity v is related to the momentum by v = 
p/E [13j . In the non-relativistic regime, when E m m 
and p ~ mv, the density p of the dark matter is given by 
p = mn, where n is the particle number density, while 
its pressure P can be obtained as [l3T | 



P 



giving 



.9 

3h 3 



1 2 
p c 

~E~ 



,f(p)d 3 p w Air 



P = v 2 P, 



3h 3 



p- 



dp, (18) 



(19) 



where a 2 = (v 2 )/3, and (v 2 ) is the average squared ve- 
locity of the particle, a is the one-dimensional velocity 
dispersion. In the non-relativistic approximation given 
by Eqs. (jig)) and ([T^| . the velocity dispersion a is a 
constant only for the case of the non-degenerate ideal 
Maxwell-Boltzmann gas, whose Newtonian analogue is 
the isothermal sphere. In the following we will con- 
sider only the specific case of the ideal gases in the non- 
relativistic regime, and we will explicitly assume that the 
dispersion velocity is a constant. 

We model the galactic dark matter as a mixture of 
two fluids, with energy densities pi,i = 1,2, constant 
velocity dispersions a^i = 1,2, pressures Pi,i = 1,2, 
and with four- velocities and W^ 1 , respectively. Since 
the two fluids have different four-velocities, we may write 
U^W^ = 1 + 6/2, where generally b is an arbitrary func- 
tion of the radial coordinate. The functional form of b 
can be obtained from Eq. (j4j), which gives 



U^W„ = - tan 2a 



Pi +Pi 
Pi +P2 



P2 +P2 



Pi +Pi 



while b can be obtained as 



tan 2a 



Pi +Pi 
P2 +P2 



P2 +P2 



Pi +Pl 



2. 



(20) 



(21) 



In the following we choose our coordinates such that 
we are comoving with the fluid with velocity U 11 . Thus 



4 



U° = e~"/ 2 and U 1 = U 2 = U 3 = 0. We also assume 
W 2 = W 3 = 0. From the conditions U^U^ = 1 and 
W^W = 1 we obtain W° = (1 + 6/4)e^/ 2 and W 1 = 
y/(b/2)(l + b/8)e~ x / 2 . By assuming that the physical 
parameters of the two fluids satisfy the condition 



2b 1 



(pi + Pi)(p 2 +P 2 ) 
8 J (Pi+ P2 + P1+P2) 2 



« 1, 



(22) 



the energy density, and the radial and tangential pres- 
sures can be obtained as 

1 , bf b\ (pi + Pi)(p 2 + P 2 ) 

£ = Pl+/'2 + -|l + - 7 j : : r, 23 

2 V 8 / KPi + P2+P1 +P2) 



b ( b 

^ = Pl+P2 + ^[l + - 



(Pl + Pi) ( P 2 + P2) 

8 J (P1+P2+P1+P2)' 



U = P 1 +P 2 



(24) 
(25) 



In the following we assume that the masses mi and m 2 
and the number densities n\ and n% of the dark mat- 
ter particles satisfy the scaling relation m\lm% = n 2 /ni, 
giving pi = p 2 = p for the densities of the two compo- 
nents. Taking into account that Pi — cr 2 pi, P2 = c r |/ ? 2 ) 
the expression of b can be written as 



b = tan 2a 



- 2. 



(26) 



Hence for this particular case it follows that b — constant. 
Moreover, since the rotation angle is arbitrary, for ai > 
U2, we can always find an angle a so that the condition 
b > is also satisfied. In the following we will restrict 
our analysis to the range of positive values of b. 

Therefore the energy density and the pressures of the 
two fluid dark matter mixture can be written as e = eop, 
* = * p, n = n p, where e = [2 + (6/2) (1 + 6/8)], 
*o = [<j\ + <j\ + (6/2) (1 + 6/8)] , and n = a 2 + a 2 , 
respectively. By introducing a set of dimensionless vari- 
ables (6,t;,m), defined as e = s c 8, r — ^/ v / 47re c , M = 
mf^/AirSc, where e c is the central density of the dark mat- 
ter halo, the continuity equation, the TOV equation and 
the tangential velocity can be written in a dimensionless 
form as 

dm (27) 



(1 



le)0(e0/le 



2m/0 

_ eo/ie 



2 (7n - 1) 



J tg 



£(l-2m/0' 



(28) 



(29) 



where j £ = eo/^o and 7n = Hq/^q, respectively. Equa- 
tions (|2T|) and (|2"8")l must be integrated with the boundary 
conditions 9 (0) = 1 and m(0) = 1. The variation of the 
two-fluid dark matter density profile and of the tangen- 
tial velocity of test bodies is represented, for fixed values 
of a 1 and 02, in Figs. [1] and [5] 
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FIG. 1: Dimensionless density profile of a two-component 
dark matter fluid, with o\ = 950 km/s and 02 = 750 km/s, 
for different values of b: b = 0.0001 (solid curve), b = 0.0002 
(dotted curve), b = 0.0004 (dashed curve), and 6 = 0.0005 
(long-dashed curve), respectively. 
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FIG. 2: Tangential velocities of test bodies in a two- 
component dark matter fluid, with <ti = 950 km/s and 
(72 = 750 km/s, for different values of 6: b — 0.0001 (solid 
curve), b = 0.0002 (dotted curve), b = 0.0004 (dashed curve), 
and b = 0.0005 (long-dashed curve), respectively. 



V. COSMOLOGICAL IMPLICATIONS 

Since the two-fluid dark matter mixture represents 
an anisotropic fluid, on a cosmological scale the corre- 
sponding geometry is also anisotropic. In a homogenous 
universe a two-fluid mixture can be described by a flat 
Bianchi type I geometry, with the line element given by 



ds 2 = dt 2 - a\{t)dx 2 - a 2 2 {i)dy 2 - a\(t)dz 2 



(30) 



where a,i(t), i = 1, 2, 3 are the directional scale factors. In 
this geometry the gravitational field equations take the 
form 



3H + Hi + Hl + Hi 



1 



(e + * + 2n) , (31) 
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Id, 



2^ 



1—(VH 3 ) = -(e-IS) 



(32) 
(33) 



where we have denoted V = 010203 , Hi = hija^i = 

1,2,3 and H = \ J2i=i H i = ^/ 3 ^> respectively. On 
a cosmological scale the pressures of the two-fluid dark 
matter obey the conditions ^,11 <C £, and hence in the 
gravitational field equations we can neglect \& and II with 
respect to e. By adding Eqs. (|3"2")l we obtain 



H + 3H 



1 



(34) 



leading to Hi = H+Ki/V, i = 1,2, 3, where K i: i — 1,2,3 
are constants of integration satisfying the consistency 
condition £ 3 =1 *Q = 0. We also denote K 2 = £ 3 =1 K 2 . 
From the conservation of the energy density of the dark 
matter, e + 3He = 0, we obtain e — ec/V, where ec 
is an arbitrary integration constant. Eq. ([3~4")l takes the 
form V = (3/2)ec, with the general solution given by 
V(t) = {i/A)e c t 2 + C x t + C 2 , with d,i = 1,2 arbi- 
trary integration constant. The large-time behavior of 
an anisotropic cosmological model can be obtained from 
the study of the anisotropy parameter A, defined by 



A = 



1 3 

i=l 



Hj-H 
H 



m 2 
v 2 



(35) 



If A = the cosmological model is isotropic. For the 
two-fluid dark matter model A = 3K 2 / [(3/2)e c * + Ci] 2 , 
and therefore in the large-time limit A 0. Thus, on 
a cosmological scale, and in the large time limit, the ini- 
tially anisotropic dark matter fluid mixture behaves like 
an isotropic fluid. 



VI. DISCUSSIONS AND FINAL REMARKS 



to two non-interacting pressureless fluids. For this case 
the theoretical model reduces to a single anisotropic fluid 
with vanishing tangential pressure. In the general case 
of arbitrary fluids there is a general r-dependent func- 
tional relationship between the energy density and the 
radial and tangential pressures, which differs radically 
from the simple barotropic equation of state. On the 
other hand, in was shown in [8] that half of the 14 galax- 
ies considered in the study are well fit by the polytropic 
halo model, despite its serious physical simplifications. 
The necessity of considering non-standard dark matter 
models with pressure is justified by our uncertainties in 
the knowledge about the nature of the dark matter par- 
ticle, as well as by the fact that this model gives a much 
better description of the observational results, as com- 
pared to the pressureless case. As a consistency check of 
our model we have considered the cosmological behavior 
of a mixture of two fluids, and we have showed that such 
a mixtures always isotropizes in the large time limit. 

The approach used in the present paper is different 
from the method used in [ll[ for the study of dust shells. 
First of all, the two fluids are described as a single ef- 
fective anisotropic fluid. While in [ll[ a non-comoving 
frame of reference is adopted, due to our transformation 
of the energy-momentum tensor, we can write down and 
study the field equations in a comoving frame. More- 
over, the inclusion of the pressure terms is straightfor- 
ward. We have also to point out the differences in geom- 
etry: while in fll| the time-dependent Lemaitre-Tolman 
model is used, and additional symmetries are imposed in 
order to reduce the number of variables, in the present 
paper we are considering the properties of dark matter 
in a spherically symmetric static geometry. The present 
model can be easily generalized for other equations of 
state of dark matter, or for cases involving the presence 
of condensates or scalar fields. These possibilities, and 
their impact on dark matter distribution, as well as some 
cosmological implications, will be considered in some fu- 
ture studies. 



In the present paper we have considered the theoretical 
possibility that dark matter may be modeled as a mix- 
ture of two non-interacting perfect fluids, with different 
four-velocities. This model is formally equivalent to a 
single anisotropic fluid. By considering a purely kinetic 
description of dark matter the basic equations describ- 
ing the dark matter mixture have been obtained, and 
integrated numerically. In the particular case in which 
$,!!<£, we can approximate "J = II « 0, corresponding 
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